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Abstract. We apply the functional renormalization group method to the calculation 
of dynamical properties of zcro-dimcnsional interacting quantum systems. As case 
studies we discuss the anharmonic oscillator and the single impurity Anderson model. 
We truncate the hierarchy of flow equations such that the results are at least correct 
up to second order perturbation theory in the coupling. For the anharmonic oscillator 
energies and spectra obtained within two different functional renormalization group 
schemes are compared to numerically exact results, perturbation theory, and the 
mean field approximation. Even at large coupling the results obtained using the 
functional renormalization group agree quite well with the numerical exact solution. 
The better of the two schemes is used to calculate spectra of the single impurity 
Anderson model, which then are compared to the results of perturbation theory and the 
numerical renormalization group. For small to intermediate couplings the functional 
renormalization group gives results which are close to the ones obtained using the 
very accurate numerical renormalization group method. In particulare the low-energy 
scale (Kondo temperature) extracted from the functional renormalization group results 
shows the expected behavior. 
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1. Introduction 

The reliable calculation of physical properties of interacting quantum mechanical 
systems presents a formidable task. Typically, one has to cope with the interplay of 
different energy-scales possibly covering several orders of magnitude even for simple 
situations. Approximate tools like perturbation theory, but even numerically exact 
techniques can usually handle only a restricted window of energy scales and are 
furthermore limited in their applicability by the approximations involved or the 
computational resources available. In addition due to the divergence of certain classes 
of Feynman diagrams some of the interesting many-particle problems cannot be tackled 
by straightforward perturbation theory. 

A general scheme that is designed to handle such multitude of energy scales and 
competition of divergences is the renormalization group pQ. The idea of this approach is 
to start from high energy scales, leaving out possible infrared divergences and work one's 
way down to the desired low-energy region in a systematic way. The precise definition 
of "systematic way" does in general depend on the problem studied. 

In particular for interacting quantum mechanical many-particle systems, two 
different schemes attempting a unique, problem independent prescription have emerged 
during the past decade. One is Wegner's flow equation technique 0, where a given 
Hamiltonian is diagonalized by continuous unitary transformation. From the final result 
one can extract detailed information about the structure of the ground state and low- 
lying excitations. This technique has been applied successfully to both fermionic and 
bosonic problems [3] . However, in general it becomes a rather cumbersome task to really 
calculate physical quantities, especially dynamics from correlation functions. Here, one 
typically has to introduce further approximations [UE|, again tightly tailored for the 
problem under investigation. 

The second field theoretical approach, which we want to focus on in the following, 
is based on a functional representation of the partition function of the system under 
consideration. It has become known as functional renormalization group (fRG) 
HI El M ■ A detailed description of the fRG will be given in the next section; here we 
make some principle remarks and discuss previous applications. 

Different versions of the fRG have been developed over the last few years. One either 
generates an exact infinite hierarchy of coupled differential equations for the amputated 
connected m-particle Green functions of the many-body system O IH] , or the one-particle 
irreducible m-particle vertex functions [3 |H] respectively. For explicit calculations one 
has to truncate the set of equations which is the major approximation involved. At 
what level this truncation is performed to obtain a tractable set of equations depends 
on the complexity of the problem. We here exclusively study the one-particle irreducible 
version of the fRG. It has the advantage of including self-energy corrections already in 
low truncation orders and being formally easily extendable to higher orders jlOj . 

Up to now, most applications of the fRG in many-body physics concentrate on 
low-dimensional, interacting fermion systems where it provides a possibility to sum 
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diverging classes of diagrams. The homogeneous two-dimensional Hubbard model 
ITB*] . the homogeneous one-dimensional Tomonaga-Luttinger model [Tl], and 
the one-dimensional lattice model of spinless fermions with nearest neighbor interaction 
and local impurities jT^I were investigated. The focus was put on properties of the 
system close to the Fermi surface, for example on the hierarchy of interactions to identify 
possible instabilities [TTJ [TJJ HB] or on Tomonaga-Luttinger liquid exponents [TH ITS']. 
The frequency dependence of the vertex functions was mostly neglected |16j . 

In this paper we investigate the frequency dependence of the self-energy and the 
effective interaction. For this purpose, we study two different zero- dimensional (local) 
models: the quantum anharmonic oscillator and a well-known problem of solid state 
physics, the single impurity Anderson model (SIAM). The former has quite often been 
used as a "toy model" to investigate the performance of different approximation schemes 
of many-particle physics E| • For this problem conventional perturbation theory is 
regular — although it generates a generic example of an asymptotic series [IT] - and 
one expects that compared to perturbation theory the fRG leads to a better agreement 
with the exact solution at larger coupling constants ("renormalization group enhanced 
perturbation theory"). We calculate the ground state energy and the energy of the first 
excited state as well as the spectral function of the propagator. Exact results for these 
observables can numerically be obtained quite easily. 

The SIAM has a known hierarchy of energy scales, and presents a true challenge to 
any many-body tool due to the generation of an exponentially small energy scale, the 
Kondo scale, leading for example to a sharp resonance in the single-particle spectrum. 
No exact solutions for dynamical quantities of the model are known. We present 
fRG results for the one-particle spectral function of the impurity site and compare 
them to conventional second order perturbation theory (in the interaction U of the 
impurity site) and Wilson's numerical renormalization group (NRG). In both models 
the truncated fRG scheme, which is correct at least to second order perturbation theory 
in the interaction, leads to a considerable improvement compared to plain second order 
perturbation theory. 

The paper is organized as follows. In the next section we present a detailed 
discussion of the fRG. The third part contains the application to the quantum 
anharmonic oscillator, while in the forth section we discuss the fRG for the single 
impurity Anderson model. A summary and outlook concludes the paper. 

2. Functional renormalization group 

Expressed as a functional integral the grand canonical partition function of a system of 
quantum mechanical particles (either fermions or bosons) interacting via a two-particle 
potential can be written as jTH] 




(1) 
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with 



(2) 



and Z being the non-interacting partition function. Here i\> and ip denote either 
Grassmann (fermions) or complex (bosons) fields. The multi-indices kj stand for the 
quantum numbers of the one-particle basis in which the problem is considered (e.g. 
momenta and spin directions) and Matsubara frequencies u. We have introduced the 
short hand notation 



k,k' 



with the propagator Q° of the related non-interacting problem given as a matrix. 
The anti-symmetrized (fermions) or symmetrized (bosons) matrix elements of the two- 
particle interaction are denoted by v k ' iik ^ klik2 . They contain the energy conserving factor 
o~u+u',v+u' and the factor 1/(3, with (3 — 1/T being the inverse temperature. We consider 
units such that h = 1 and ks — 1- The generating functional of the m-particle Green 
function is given by 



1 



-l 



-fan)-^)} , (3) 

with (i/i, r^j = J2k ^kVk and the external source fields n and fj. From this the generating 
functional of the connected m-particle Green function follows as 

W c (M,M)=ln[W({r7},M)]. (4) 

The (connected) m-particle Green function G$ can be obtained by taking functional 
derivatives 

^ro (X' • • • ' k'm'i 

By a Legendre transformation 



■ i km) 



$Vk[ ■ ■ ■ 8r} k > m 8n km . . . Sn kl 



W^({fj},{r,}) 



.(5) 



r)=0=ri 



r ({0}, {0}) = -W c ({f}}, {v}) - (0, V) - (V, 0) + 0, [G°] " 



(6) 



the generating functional of the one-particle irreducible vertex functions 7 m , with 
external source fields 6 and and 



7m (^1; • • • ) ^mi ^1) 



i km) 



(? 



r {0},{0} 



(7) 



<f>=0=<j> 



is obtained. Note that in contrast to the usual definition JH] of V we have added a term 
(0, [G ]' 1 ^) in Eq. © for convenience (see below). The relation between the G$ and 
7 m can be found in text books [TP] . The O-particle vertex 7 provides the interacting 
part of the grand canonical potential Q 



n 



-TlnZ = T70 — T In Zq 
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For the 1-particle Green function we obtain 
G 1 (k' ] k)=G c 1 (k';k) = -(g k ,, k 

where 



7i - C 



-i 



-i 

k'k 



Ik'k 



— E 



k'k 



with the self-energy £, and ( = — 1 for fermions or ( = I for bosons, respectively. 
This implies the relation £ = £71. The 2-particle vertex 72 is usually referred to as 
the effective interaction. For m > 3, the m-particle interaction 7 m has diagrammatic 
contributions starting at m-th order in the two-particle interaction. 

In Eqs. (JIJ and Q we replace the non-interacting propagator by a propagator Q 0,A 
depending on a cutoff A and Z$ by Z A determined using Q 0,A . The boundary conditions 
for the cutoff A £ [A , 0] are taken as 

£0,A 0=0 ) £0,A=0 = £0 ; (g) 

i.e. at the starting point A = Ao no degrees of freedom are "turned on" while at A = the 
cutoff independent problem is recovered. In our applications we use a sharp Matsubara 
frequency cutoff with 

0°- A = 6 (|u;| - A) 0° (9) 

and consider A — > 00 |T^]. Through Q 0,A the quantities defined in Eqs. (JH) to (JTj) acquire 
a A-dependence. One now derives a functional differential equation for T A . From this, 
by expanding in powers of the external sources, a set of coupled differential equations 
for the 7 A is obtained. 

As a first step we differentiate W C ' A with respect to A, which after straightforward 
algebra leads to 



d_ 

dA 



W 



c.A 



C Tr (SV' A ) +Tr ( Q 



5fj5n 



+ c 



'SW 



c.A 



5rj 



Q 



6W 



5'ij 



with 



Q: 



dA 



■>0,A 



(10) 



'11' 



Considering and <fi as the fundamental variables we obtain from Eq. © 



dA 



a' 



Applying the chain rule and using Eq. (fTUjl this leads to 



>Ar>0,A 



Tr Q 



5 2 W 



c,A ' 



5n A 5rf 



The last term in Eq. (jUJ) is exactly cancelled, which a posteriori justifies its inclusion. 
Using the well known relation ^H] between the second functional derivatives of V and 
W c we obtain the functional differential equation 
d 



dA 



r A = -c Tr (g A £°' A ) - Tr (s A v^(r A , g°' A ] 



(12) 
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where Vl' 1 , stand for the upper left block of the matrix 



A r>A\ 



V 



3 2 r A 

<50<50 



s 2 r A 



-i 



<5 2 r A 



?0,A 



(13) 



and the upper index t denotes the transposed matrix. To obtain differential equations 
for the 7^ which include self-energy corrections we express Vg <* in terms of Q A instead 
of Q A '°. This is achieved by defining 

5 2 T A a 



7i 



and using 



g 



0,A 



C7f 



i -i 



which leads to 



dA 



-( Tr (QV' A ) +CTr (Q A Q A V/AV A , £ 



'A/->A-r,l,l/-pA r-A^ 



with 





<5 2 r^ 



04 



-1 



(14) 



(15) 



(16) 



* 2 r A and g£ are at 



For later applications it is important to note that ^ as well as 
least quadratic in the external sources. The initial condition for the exact functional 
differential equation ()15j) can either be obtained by lengthy but straightforward algebra 
not presented here, or by the following simple argument: at A = A , Q°' A ° = (no 
degrees of freedom are "turned on" ) and in a perturbative expansion of the 7^° the only 
term which does not vanish is the bare 2-particle vertex. We thus find 

r Ao ({0},{0})=5 int ({0},{0}) . (17) 

An exact infinite hierarchy of flow equations for the 7^ can be obtained by 
expanding Eq. ()16p in a geometric series and T A in the external sources 

00 Am 

rA (W.W) = Er^ E E 7 m (^,...,C^---,U 



in l m 2 k\ />' 5, /.-, 



X 0« 



The equation for 7^ reads 

^7o A = -C Tr (QV' A ) + C Tr (Q A g A 

Via Q A the derivative of 7^ couples to the one-particle self-energy. For the flow of the 
self-energy we obtain 



4rT}W k) = C^r^ = Tr (S A 7 2 A (fc' ? 



dA liy ' ' ' dA 
with the so-called single scale propagator 

s A = g A Q A g A . 



(19) 
(20) 



Functional renormalization group approach to zero- dimensional interacting systems 7 

Here 7 A {k' , . . . ; k, . . .) is a matrix in the variables not explicitly written, i.e. 

72 (k', q'; k, q). Diagrammatically Eq. (fT^j) is shown in Fig. ^ 



7 2 A (fc' 

The derivative of 7^ is determined by 7^ and the 2-particle vertex 7^. Thus an equation 
for 7^ is required 
d 



\ k. 



9'. 9 

„A : 



^-7 2 A (^, k u k 2 ) = Tr (S A ^(k[, k' 2 , . . . ; fc ls fc 2) 



+ CTr(S A 72 A (...,...;fci,fc 2 ) 
+ C Tr (<S A 7 2 A (^, . . . ; h, . . .)g A 1 A {k' 2) ■ ■ ■ ■ • •) 
+ C [k[ <-> fc£] + C [fci <-> fe] + [A;i «-> 4, fei ^ fe] ) ■ (21) 

The corresponding diagrammatic representation is shown in Fig. El The right hand side 
(rhs) of Eq. (J2*T|) contains 7 A , 7 A , and the 3-particle vertex 73. For m > 1 the equation 
for ^7^ contains 7 A / with m' = 1, 2, . . . , m, m + 1. The initial condition for the 7 A ° 
can be obtained from Eq. (j!7|) and is given by 

^°(K, K\ h, k 2 ) = v KAMM , = for m ^ 2 . (22) 
We here refrain from explicitly presenting equations for ^7^ with m > 3 since later 



G A '^{k^k'r 



21 



•)) 



on the set of differential equations is truncated by setting 73 = 73 
that vertices with m > 3 do not contribute. 



An 



0, which implies 



w 

Figure 1. Diagrammatic form of the flow equation for 7^ = ££ A . The slashed line 
stands for the single scale propagator S A . 



^2 * ^2 ^2 ^""^^ ^2 ^2 

ki k-i ki k-i k-i 



k? 



ki 




ki wk. 




kzw^ 



^Aki k 2 Ak2 k,Ak2 k 2 Ak; 



Figure 2. Diagrammatic form of the flow equation for 7^. The slashed line stands 
for the single scale propagator S A the unslashed line for Q A . 



Following the above systematics, a truncation scheme emerges quite naturally. If, 
for m c > 2, the vertex 7 A c+1 on the rhs of the coupled flow equations is replaced by its 
initial condition 7 m ° , x = 0, a closed set of equations for 7^ with m < m c follows. This 
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set of differential equations can then be integrated over A starting at A = A down to 
A = providing approximate expressions for the 7 m of the original (cutoff free) problem 
with m < m c . Expanding 7^ in terms of the bare interaction, conventional perturbation 
theory for the grand canonical potential, the self-energy, the effective interaction and 
higher order vertex functions can be recovered from an iterative treatment of the flow 
equations. It is easy to show that the vertex functions obtained from the truncated 
equations are at least correct up to order m c in the bare interaction. 

To obtain a manageable set of equations in applications of the fRG to one- and two- 
dimensional quantum many-body problems fTJ EIJ El E] > further approximations 
in addition to the above truncation scheme were necessary. In the following two sections 
we will avoid such additional approximations and solve the truncated fRG equations to 
order m c = 2 for two zero-dimensional (local) interacting many-particle problems: the 
quantum harmonic oscillator with a quartic perturbation and the SIAM. 

Very recently a modified version of the flow equation for 7^ was suggested. As we 
will consider the truncation order m c = 2 we here only describe this modified scheme 
for this order. Guided by the idea of an improved fulfillment (compared to the scheme 
described above) of a certain Ward identity Katanin replaced the combined propagator 
S A in the last five terms of Eq. (|21|) (the first term does not contribute to order m c = 2) 
by -dQ A /dk [H|. Using 

dQ A * a dYj A a 

— — — = s A + q a — q a 

dk dk 

obtained from Eq. (|Tl)l it is obvious that the terms added to Eq. (|2*T|) are at least of third 
order in the bare interaction. Using this modification in the order m c = 2 equations 
for 7^ and 7^ leads to the exact solution for certain exactly solvable models (reduced 
BCS model, interacting fermions with forward scattering only). This suggests that the 
replacement possibly improves the results of the truncated fRG also for other models. 
In the next section we show that this is indeed the case for the harmonic oscillator with 
a quartic perturbation. The same holds for the SIAM as discussed in Sect. 4. 



3. Application to the anharmonic oscillator 

In appropriate units the Hamiltonian of the harmonic oscillator with a quartic 
perturbation is given by 

H = \x 2 + \ P 2 + fx 4 , (23) 

with the position operator x, the momentum operator p, and the coupling constant g. 
We here focus on T = and are interested in low-lying eigenenergies E n as well as the 
(imaginary) time-ordered propagator 

g(r) = (E \T[x(r)x(0)}\E ) , 
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with \E n ) being the eigenstates of the Hamiltonian in Eq. (|23|). The Fourier transform 
of the propagator can be written as 

1 



Q{iu) 



dre luJT G( 



T 



where we have introduced the self-energy £ and the non-interacting propagator 



LJ 2 + 1 



In contrast to the more general notation used in the last section, propagators and the 
self-energy only depend on a single frequency and do no longer contain the energy 
conserving ^-function (T = 0) here. The propagator has the Lehmann representation 



Qiiuj) 



-| oo 

-Y 

n=l 



1 



1 



(E \(a + a^\E n ) 



iu + (E n - E Q ) iu - (E n - Eq 
where a, a) denote the usual raising and lowering operators. The spectral weights 



(24) 



(25) 



(E \(a + a^ \E n ) 
and energies E n fulfill the f-sum rule 

oo 

l = Y i (E n -E ) s n . 

n=l 

It turns out that for coupling constants g < 50 considered here the sums in Eqs. (|24|) and 
(I25|) are dominated by the first term. For this reason only the first few eigenstates and 
eigenenergies are required to obtain accurate ("numerically exact") results for Q{iuj). 
These can quite easily be obtained by expressing H in the basis of eigenstates |n) of the 
unperturbed harmonic oscillator and numerically diagonalizing the upper left corner of 
the (infinite) matrix (n \ H \ n') with n, n' < n c and a sufficiently large n c . For g < 50, 
n c = 100 turns out to be large enough to fulfill the sum rule Eq. (|23jl to very high 
precision. 

Second order perturbation theory for the ^-dependent part of the ground state 
energy yields 



,( 2 ) 



E, 



(2) 



1 



7 



32 

and for the self-energy one obtains 



1536 



9 



1 1,1 

- q H Q + - 

4 y 32 y 8 



9 



(26) 



(27) 



Within the fRG approximate expressions for 



E n n — E n — E 



and s n can be obtained from the poles and residues of the propagator Q{iuj). 
Furthermore, since Eqs. ()24j) and (|2"K|) are dominated by the first terms we only consider 
Eifi and s± [21]. Second order approximations for these quantities are given by the 



smallest pole of Q^'{ 



+ 1-E( 2 )( 



ioj , 



and the related residue. It is important 



(2) (2) 

to note that this approximation E{ q agrees with E{ 



E^\ where E± is determined 
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directly from Rayleigh-Schrodinger perturbation theory, only up to second order in g, 
but is closer to the exact E li0 . 

Within mean field theory, e^ F and a frequency independent E MF are given by 

er = - f <* 2 >l - \ ■ sMF = -f (^ 2 >m, 

and (X 2 ) MF is the solution of the self-consistency equation 



I MF 



<l + g(X*) MF /2 
From the mean field propagator one obtains 

£f F =Jl- S MF S MF = I^^ = _ 

i,o v ' 1 2 a/1 - S MF 

The functional integral representation of the grand canonical partition function of 
the Hamiltonian Eq. (}2*3*|) reads 
Z 1 



— = —JVxxexp! [ (x,[g } 1 x)/2-S mt ({x},{x})} , (2* 



with 



Sint {{x}, {x}) = — ^ni+n 2 +n 3 +n 4 ,0 x{iUx)x{iu 2 )x{iu^)x{iu i ) , (29) 

ni ,...,714 

bosonic Matsubara frequencies ox,- = 27inj/j3, and complex fields ic(ia>) = x(— iw). As 
outlined in the last section and using a frequency cutoff Eq. Q flow equations for the 
7 A can be obtained. Here we focus on the equations in truncation order m c = 2. For 
T — > and after introducing 



we find 



e A = lim T 7o A 



d A 1 , 

-rren = m 

dA 2vr 



£°(zA)£ A (A) , (30) 



with the initial condition e A= °° = 0. At the end of the flow, e A=0 directly provides the 
fRG approximation eo RG for the (^-dependent part of the ground state energy. The flow 
equation for the self-energy follows as 

^r sA M = t- aui 1 vAr S ^ -iu>, iA, -iA) , (31) 
aA 2n A z + 1 — L A (zco>) 

with the initial condition £ A= °° = and the fRG approximation for the self-energy 
£ fRG (ia;) = T, A=0 (iuj). Here g A denotes the totally symmetric 2-particle vertex which, 
in contrast to the vertex 7 A introduced in the last section, does not contain an energy 
conserving 5-function and factors of (3. It depends on only three frequencies, but the 
fourth will nevertheless always be included in the following. To derive Eqs. (J3*U|) and 
(131(1 one has to deal with products of delta functions 5(\u\ — A) and terms involving 
step functions 0(|u;| — A). These seemingly ambiguous expressions are well defined and 
unique if the sharp cutoff is implemented as a limit of increasingly sharp broadened 
cutoff functions e , with the broadening parameter e tending to zero. The expressions 
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can then be conveniently evaluated by using the following relation [8 , valid for arbitrary 
continuous functions /: 

6 e (x - A) f[G e (x - A)] -> 8(x - A) f 1 f(t) dt , (32) 

where 5 t = —dQ e /de. 

For g A the flow equation reads 



d I f°° 

—g A (iu lj iu 2j i&3, —uo\ — iu 2 — i^s) = tt~ / d v [P(iv, iv — iu\ — iu 2 ) 
dA zir J-oo 

~A/- . • . • ... • . • . \ A 



x g (iu>i,iu>2i —iv, iv — iu\ — iu%)g (iu>3, —iu\ — iu 2 — iu>s, —iv + iu% + iuj^-, iv) 
+ (lu 2 <-> uj 3 ) + (lu 2 <-> — UJ\ — uj 2 — Us)] (33) 

where V(iv,ii/) stands for two different products of propagators. In the conventional 
fRG scheme it is given by 

V cm {iv^v')=S K {iv)Q K {iv') (34) 

while in the modified scheme j2D] one obtains 

V mod (iv, iv') = -^P- £ V) • (35) 
dA 

To explicitly evaluate V(iv, iv') Eq. has to be used, where special care has to be 
taken for the case v = v'. In the conventional scheme V CO n{iv,iv') contains a factor 
S(\v\ — A) and the integral over v in Eq. (J3*3*j) can be performed analytically. 

To numerically solve the set of differential equations (jSOJ), dSJ), and (|3~3|) we have 
discretized the frequencies (which at T = are continuous) on a linear mesh Uj = jS 
with j = —jo, — jo + 1, ■ ■ ■ ,jo f22j - By increasing j and decreasing 5 convergence can be 
achieved up to the required accuracy. For our purposes jo = 40 and 5 = 0.5 turned out 
to be appropriate. This leads to a set of roughly 5.3 x 10 5 coupled equations. The A- 
integration is started at Ao = 10 5 making sure that further increasing Ao does not lead to 
significant changes in the results. Figs. El to El show comparisons of e , ^.q, and sj with 
the different approximations considered here (second order perturbation theory, mean 
field, conventional fRG, modified fRG) and the exact results. Although the approximate 
fRGs correctly reproduce only the first two derivatives with respect to g at g — 0, they 
give extremely accurate results even up to g = 50, while conventional second order 
perturbation theory can only be trusted for g < 1. This provides an impressive example 
of the power of "renormalization group enhanced perturbation theory" . Comparing the 
two fRG approximations the modified scheme is roughly a factor of two closer to the 
exact result and thus indeed a substantial improvement [201 E3] ■ 

To avoid the problem of analytic continuation (see the next section) the results for 
ef^ G and s^- G were obtained by fitting a function a/(u 2 + b 2 ) with a and b as fitting 



parameters to Q (iu) = u +1 — S (iu) . Assuming this fitting form we have 
used that the spectral function is dominated by the first peak [2IJ. For the problem 
studied also mean field theory leads to fairly accurate results (but not as good as the 
fRG). This is related to the fact that low-lying eigenstates of the Hamiltonian in Eq. 
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Figure 3. Coupling constant dependent part of the groundstate energy eo = Eq — Eq 
as a function of g. Different approximations (second order perturbation theory [thin 
dashed line], mean field [dotted line], conventional fRG [thick dashed line], modified 
fRG [dashed-dotted line]) are compared to the exact result (solid line). The inset shows 
the difference between the exact result and the two fRG approximations. 



2.5 1 ■ 1 ' 1 ' 1 ■ r 




Figure 4. As in Fig.|21 but for the energy difference £1,0 of the first excited state and 
the ground state. 
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0.5 1 n 1 ■ 1 ■ 1 ■ 1 ■ 1 




Figure 5. As in Fig. but for the spectral weight si of the first peak. 

(I23|) can be described quite well by the eigenstates of a harmonic oscillator with a shifted 
frequency determined self-consistently. 

4. Single-particle dynamics of the single impurity Anderson model 

In contrast to the anharmonic oscillator studied in the previous section, the single 
impurity Anderson model [23] 

H = H es4c fc + e d £ S a d a + Udfad^ + X £ (ct^ + h.c.) (36) 

fccr CT fccr 

consists of two subsystems, namely a "conduction band" with continuous energy 
spectrum described by the first term in Eq. (J36j) and a localized level (typically referred 
to as "oT -state) with energy relative to the chemical potential \i of the conduction 
electrons. Two electrons occupying the localized level are in addition subject to a 
Coulomb repulsion U > 0. Both subsystems are coupled via the hybridization in the 
last term. As the two-body interaction is restricted to a single level the SIAM falls into 
the of class of zero-dimensional interacting systems. 

While this model looks rather simple, it contains all ingredients that make it a 
complicated many-body problem. The bare energy scales of the model are the bandwidth 
W of the band electrons, the local energy e^, the Coulomb repulsion U, and the bare 
level width generated by the hybridization, A = nV 2 p c (n), where p c (e) denotes the 
density of states of the conduction states which is assumed to be slowly varying. In the 
following we restrict ourselves to the particle-hole symmetric case 2e^ + U = 2/x = 0. 
As in the parameter regime U/Aq 3> 1 the charge fluctuations on the d-leve\ from the 
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average value 1 are small, it can effectively be described by a spin antiferromagnetically 
coupled to the conduction electron spin density [23], the so-called Kondo model [26J. 
This antiferromagnetic coupling leads to a screening of the local spin by the conduction 
electrons with a characteristic energy scale ln(T K ) oc — 1/A , the Kondo temperature. 
As is apparent from this expression, Tk depends non-analytically on Ao, signalling the 
occurrence of infrared divergences in perturbation theories in A /U j^Z] and a severe 
problem for computational techniques, namely the task to resolve an exponentially small 
energy scale. 

On the other hand, the physics in the regimes T, u 3> T K and T, uj <C T K is 
comparatively simple. For T, u ^> Tk it is governed by charge excitations with energies 
and e^ + U with a life-time given by 1/(2A ). In the other limit T,u <C Tk it has 
been worked out using Wilson's NRG [28] that the system can again be described by a 
Hamiltonian of the form Eq. (jSEJ)? but with U — > U* = 0, td — > e* d = and A replaced 
by Aq ~ Tk- This regime has been coined "local Fermi liquid" by Nozieres This 
effective description yields a narrow resonance at the chemical potential in the spectral 
function of the ci-Green function which in the symmetric case takes the form 

G*(i") = -■ xA vcY (37) 

iv — A{tv) — L(zz/J 

Here A(iv) = V 2 J dep c (e) / (iv — e) and £ denotes all self-energy contributions of 
second and higher order in U. As a special feature of the symmetric case, already the 
approximation to only keep the second order self-energy X^ 2 ) describes the qualitative 
behavior of the <i-spectral function for different values of U/A correctly [21] • For 
U/Aq <C 1 there is a Lorentzian peak at the chemical potential with a width differing 
little from A . For U/A 3> 1 most of the spectral weight is in the high energy peaks 
near ±U/2 with a narrow resonance at \l = 0. Quantitatively the width and shape of 
this Kondo (or Abrikosov-Suhl) resonance is described poorly using T,( 2 \ vanishing only 
~ Aq/U. From the results of Sect. 3 one expects that the application of the functional 
renormalization group to the SIAM leads to improvements over the direct perturbational 
description of the Kondo resonance. For the SIAM the fRG results can be compared to 
the outcome of NRG calculations PJEHJ- An important aspect of the treatment of the 
SIAM with the fRG is that, while the numerical effort of the NRG method increases 
exponentially with the number of impurity degrees of freedom, e.g. in case of additional 
orbital degrees of freedom or for a system of many magnetic impurities, the increase in 
computational resources necessary in the fRG approach described below is governed at 
most by a power-law. 

The starting point of our fRG approach are again Eqs. (fTUj) and ()21|) . We also use 
m c = 2, i.e. we replace the 3-particle-vertex by its initial condition 7^ = 0, and consider 
T = only. We again use the frequency cutoff described in Eq. Q and because we 
presently concentrate on spectral properties only flow equations for the self-energy and 
for the 2-particle-vertex (and not the ground state energy) will be derived. Note that 
the calculation of two-particle properties is possible within the same scheme without 
further difficulties ^2] and will be discussed in a forthcoming publication. 
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For the SIAM all Green functions Gij different from Gdd can be expressed in terms 
of Gdd and Green functions of the non-interacting system. For the calculation of Gdd 
the multi-index k is replaced by a spin-index a and a frequency v. The 2-particle-vertex 
can then be written as 

7 2 A (^, k' 2 ; fci, k 2 ) = Sfa + v 2 -v[- v' 2 ) (38) 
xjXi^aa,^ U A (if[,if' 2 ;ifi,iv2) -5^5^ U A (if 2 ,if[;ifi,if 2 )}. 

This can be used to perform the sum over the spins in Eq. ()19j) and we find for the 
self-energy 



a A Z7T 



2U A {if, iv'\ iu, iu') - U A {if', iu; iu, iu') . (39) 



In order to derive a flow equation for IA (iu[, ii/ 2 ; if x, if 2)1 with v 2 = f[ + u 2 — fx, the 
spin-sum in Eq. (12 1|) has to be performed as well. This leads to the lengthy expression 
Eq. (jA.l|) . presented in Appendix A. 

The initial values are given by £ A= °°(iz/) = and U A=00 (if[,if 2 ,ifi,if 2 ) = U. 
In the numerical solution of these flow equations we start integrating at finite Ao ^ 
max ({7, A ). The integration from A = oo to A can be performed analytically. 
Up to corrections of order Aq 1 the new initial conditions are T, A °(iu) = U/2 and 

U °(if' 1 ,if 2 ;ifi,if2) = U. 

S A (if) and V A n (if,if') = S A {if)Q A {if') are calculated the same way as for the 
anharmonic oscillator using Eq. (|32|) . and are given by 

5 A H - *(M " A) ] , (40) 



with [£/°(«V)] 1 = if — A(iv) — €d and 



V A n (if, if') -> <5(|H - A) ^ 9(1 ^ A) (41) 

with O(0) = 1/2. In our calculations we use the limit of an infinite bandwidth, i.e. 
A(iu) — ► — i sign(z/) A . The results of second order perturbation theory for the self- 
energy can be recovered by replacing the self-energy on the rhs of Eqs. (jBHJ) and (|A.1J) and 
U A on the rhs of Eq. (|A.1|) by their initial values. It turns out that the fRG version using 
V A n gives good results for U/A < 3, but fails to provide the expected improvement 
compared to the use of E^ 2 -* for U/Aq > 3. This will be discussed in more detail in 
a forthcoming publication. For this reason we follow [2U| and replace V^ n (if,if') by 



V A od (if,if') = -g A (if')-^G A (if) in Eq. (J23J). Applying Eq. © yields 
>A,^ d nKfi ^ , x , u>| A , 1 0(|i/| -A) 



■g A (!i/')-g A («^) -> <%|-a) 



rfA v 7 Vl 1 [^°(^)] -£ A (w) [£°(zi/)l -£ A (ii/) 

d sA e(H-A) e(H-A) 



(42) 



dA y ' [[g°{iu))- L - £ A (iz/)] 2 [G (ij/)]~ - E A {iv 
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For the integration of the T = 0-flow equations the continuous frequencies again 
have to be discretized. However, in contrast to the anharmonic oscillator it is not 
sufficient to work with a linear mesh. Because for U/Aq ^> 1 the Kondo resonance 
is visible only on an exponentially small energy scale around the Fermi level, we use 
a combination of a linear and a logarithmic mesh. This enables us to recover both, 
the high energy as well as the low-energy physics keeping the number of frequencies 
to a manageable size. The numerical effort grows with the third power of the number 
of frequencies in the conventional version and with almost the forth power using the 
modified version, due to the absence of a 5-function in the last term of Eq. (J35J). 

In the fRG one naturally obtains the self-energy for imaginary frequencies. In order 
to calculate spectral functions, an analytic continuation to the real axis is necessary, 
which is known to be an ill-posed problem. Since the results from fRG are not subject to 
statistical errors or noise as for example quantum Monte Carlo data, we have applied the 
method of Pade approximation |2H| to obtain from H{iv). We find that especially 
the low-energy part of the spectral function which contains the Kondo resonance can be 
reliably extracted if sufficiently many frequencies close to the Fermi level are used (see 
below). 




00/ A n Q)/A, 



Figure 6. Left panel: T,(ioj) from fRG for U = Ao (circles) compared to NRG (dashed 
line) and second order perturbation theory (dotted line). 

Right panel: 3mS(ti> + i0 + ) for fRG obtained from Pade approximants for the data in 
left panel. The inset shows an enlarged view around w = 0. 

As an example, we present in Fig. |B] the calculated (left panel) and the 

corresponding 3mS(u + i0 + ) from the Pade approximation (right panel) for U/Aq = 1. 
The dashed and dotted lines represent results from NRG and second order perturbation 
theory, respectively. The inset in the right panel shows an enlarged view of the region 
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around uj = 0. Apparently, the Pade approximation to the fRG provides reliable results, 
in particular around uj = 0, and recovers the perturbation theory as expected for such 
small value of U. Discrepancies to NRG for large uj can be traced to broadening effects 
in the NRG. 

The true test for the method however is a comparison of the self-energy to NRG 
results for larger values of U and in particular the behavior on low-energy scales. Such 
a comparison is presented in Fig. for U/A = 1, 5, and 10. For U/A = 5 the 
fRG results are still in excellent agreement with the NRG data while second order 
perturbation theory already deviates. For U = 10A we observe a significant dependency 
of the general structures and also of the slope at uj — > on details of the discretization 
mesh. It seems necessary to have an extremely fine resolution around uj —>■ and a 
sufficient resolution around uj ~ U/2. With an exponentially vanishing low-energy scale 
these constraints are hard to fulfill with both linear and logarithmic meshes keeping the 
number of flow equations to a manageable size. The bottom left pannel contains the 
best fRG data we were able to obtain so far. A more extended discussion of this issue 
will be presented in a forthcoming publication. 




0.1 




0.1 



-0.2 
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Figure 7. H(ioj) from fRG (circles) for U/Aq = 1, 5 and 10 compared to NRG (dashed 
line) and second order perturbation theory (dotted line). 

Bottom right: Effective mass m* from fRG (circles), NRG (squares) and second order 
perturbation theory (triangles). 
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In addition to the self-energy the effective mass m* obtained from 

* QmEjiu) 
m = 1 + hm 

as a function of U is depicted, in the bottom right panel of Fig. [7| This quantity is 
directly related to the Kondo temperature j2ZJ- Compared to perturbation theory the 
effective mass is much closer to the very accurate values determined from the NRG. 

The evolution of the spectral function of the d- level, Pd(u) = —-^mGd{u + i0 + ) 
for the above values of U is collected in Fig. |HJ The inset shows a comparison to NRG 
and perturbation theory for the region around to — at U/A = 10. One can see very 
nicely the development of the sharp resonance in the spectrum and the formation of 
the Hubbard bands with increasing U. Due to the insufficient resolution at larger oj 




Figure 8. The evolution of the fRG spectral function for the d- level for the values 
of U in Fig. The inset shows a comparison to NRG for U/Aq = 10 for the region 
around oj = 0. 

in the calculation for U/A = 10, the Hubbard bands come out too broad here. We 
believe that an improved discretization of the energy mesh in the solution of the flow 
equations will remedy that particular problem. The important region around oj = 
on the other hand is captured rather well by the fRG. It is in particular noteworthy 
that the functional form at low frequencies (see the inset) apparently does not follow a 
Lorentzian but rather, like the NRG, the more complex scaling form with logarithmic 
tails predicted by Logan et al. fTTTTj . 
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5. Summary and outlook 

We have presented an application of the functional renormalization group technique to 
solve for the dynamics of zero-dimensional interacting quantum problems. As particular 
examples, we discussed the anharmonic oscillator and the single impurity Anderson 
model. In both cases the fRG proved to be a substantial improvement over conventional 
low-order perturbation theory and rather close to the very accurate results obtained 
numerically. 

We also investigated the differences between conventional fRG and a modification 
suggested by Katanin 20J, which should improve the accuracy of the method further. 
That this is indeed true was shown directly for the anharmonic oscillator. For the single 
impurity Anderson model it was actually necessary to use this modified version to obtain 
sensible results for values U/Aq > 3. 

An important aspect of the fRG is that an extension of the SIAM to more 
complex systems, like e.g. orbital degrees of freedom or systems of coupled impurities 
is straightforward. This is in principle also true for Wilson's NRG. In the latter, 
however, the exponentially increasing Hilbert space renders a practical application 
quickly impossible. For the fRG, on the other hand, the major modification will be 
an increase of the number of equations, which means that the numerical effort will 
increase at most with a power-law. This feature makes the fRG a possible method 
to study features of complex impurity systems and in particular a potential "impurity 
solver" for mean-field theories of interacting lattice models like the Hubbard model 
in the framework of the dynamical mean- field theory 3T or the dynamical cluster 
approximation [32J. Furthermore, as it becomes clear from the anharmonic oscillator, 
the fRG is of equal complexity for bosonic and fermionic systems. This makes opens 
the possibility to study combinations of such degrees of freedom in impurity models. A 
question which can be addressed using the fRG is the influence of phonons or magnetic 
fluctuations on low energy scales. Within the dynamical mean-field theory the metal- 
insulator transition in the presence of phonons (Kondo volume collapse 33J) as well 
as the problem of non-Fermi liquid formation in for example CeC\iQ- x Au x can be 
studied. 
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Appendix A. 

Depending on the specific parameterization of the 2-particle-vertex flow equations for 
the relevant frequency dependent parts of the vertex and for the self-energy can be 
derived by using the spin conservation on the vertex when the sum over spins in Eq. 
(12 lj) is performed. For the form of the the vertex given in Eq. (jSHj) the corresponding 
flow equation for the self-energy is given by Eq. > an d the equation for the frequency 
dependent part of the vertex U A (iv[, iv 2 ; v\, iv 2 ) with v 2 — v[ + v' 2 — vi reads 




(A.l) 



—XA A {iv, iv\ + iv 2 — iv; iv\,iv 2 )lA A {iv' 2 , iv[; iv\ + iv 2 — iv, iv) 



U A (iui + iv 2 — iv, iv; iv\,iv 2 )U A {iv' x , iv 2 ; iv\ + iv 2 — iv, iv 




2U 



(iv[, iv; iv\, —iv\ + iv[ + iv)U A {iv 2 , —iv\ + iv[ + iv; iv 2 , iv) 



-U A 



(iv[, iv; ivi, —iv\ + iv[ + iv)U A {—iv\ + iv[ + iv, iv 2 ; iv 2 , iv) 



U A {iv, iv[; iv x , —iv\ + iv[ + iv)U A {iv' 2 , —iv\ + iv[ + iv; iv 2 , iv 





U A {iv, iv 2 ; iv\, —%V\ + iv' 2 + iv)U A [—ivx + iv' 2 + iv, iv[; iv 2 , iv) 



+ 




Again V stands either for P con or P mo d- 

Two other possible parameterizations of the 2-particle vertex are 



72 (& & 6i> 6) = +v 2 -v'x- ^2){<WAi,^<W; U A {iv' x ,iv' 2 ;iv x ,iv 2 ) 

+5 ai - a2 Sa[ -a' 2 (^n^W^ (ii/J, iv 2 ; iv x , iv 2 ) - 5 ai ^ a >U A {iv 2 , iv[; iv x ,iv 2 )) }. 



and 



7 2 (ti, & 6, 6) = S(ui + ^2 - v'x - v[ 



2 ){ S a ^. aim Us(iv[,iv 2 ;ivi,iv 2 ) 
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with 

Both these parameterizations lead to different set of flow equations for the self-energy 
and to two sets of flow equations for the frequency dependent functions Up and for 
the first and Ug and for the second (compared to one set for the parameterization 
Eq. (jHHjl ) implying an increased numerical effort. They are usefull to further investigate 
the processes occurring during the integration of the flow equations. This way we are 
able to distinguish the behavior of different channels of the interaction. 
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